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Control of Stimulated Raman Scattering in the Strongly Nonlinear and Kinetic 
Regime Using Spike Trains of Uneven Duration and Delay: STUD Pulses 
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Stimulated Raman scattering (SRS) in its strongly nonlinear, kinetic regime is controlled by 
a technique of deterministic, strong temporal modulation and spatial scrambling of laser speckle 
patterns, called Spike Trains of Uneven Duration and Delay (STUD pulses) [B. Afeyan and S. 
Hiiller, Phys. Rev. Lett, (submitted)]. Kinetic simulations show that use of STUD pulses may 
decrease SRS reflectivity by more than an order of magnitude over random-phase-plate (RPP) or 
induced-spatial-incoherence (ISI) beams of the same average intensity and comparable bandwidth. 

PACS numbers: 52.35.Fp, 52.35.Mw, 52.38.Bv, 52.38.Dx 



Laser-plasma instabilities (LPI) pose a risk to the real- 
ization of laser-driven inertial confinement fusion (ICF) 
ignition [1]. The present approach is to use continuous, 
ns-time-scale illumination of a target with high-intensity 
laser beams. However, because of LPI, this may prove to 
be less than ideal when compared with a novel technique 
invented by Afeyan [2] employing intermittent, scintil- 
lating, space-time illumination which may significantly 
reduce the levels of nonlinear optical processes. The ef- 
ficacy of this technique, which employs Spike Trains of 
Uneven Duration and Delay (STUD pulses), has been 
demonstrated in the fluid regime of instability evolution 
from low to moderate gains per speckle, where the lin- 
ear growth is halted by the use of STUD pulses and any 
saturation is from pump depletion |2H4j. This Letter fo- 
cuses on the application of STUD pulses to Stimulated 
Raman Scattering (SRS) in settings where kinetic nonlin- 
earity dominates the evolution of driven electron plasma 
waves (EPW) and where multi-laser-speckle, cooperative 
behavior can proceed through the exchange of hot elec- 
trons and SRS scattered light among laser speckles [5]-[7] . 
We find from our initial study that more than order-of- 
magnitude reduction in SRS reflectivity can be achieved. 
The key is to keep SRS growth below the level where 
secondary, nonlinear processes causing cooperative be- 
havior among hot spots can occur, thus disallowing the 
self-organized state. 

SRS is the resonant, three- wave coupling of a light 
wave into scattered light and electron plasma waves. It 
is an LPI process occurring in large amounts in indirect- 
drive ICF experiments with potentially deleterious ef- 
fects, including scattering of laser energy out of the 
hohlraum, redirection of energy within the hohlraum, 
and generation of hot electrons that may contribute to 
capsule preheat. At the National Ignition Facility (NIF), 
experiments show ~50% inner-cone beam energy loss 
to SRS [1 . Laser facilities such as OMEGA and the 
NIF employ beam smoothing, whereby random phase 
plates (RPP) break the laser beams into "speckles" to 
effect a quasi-uniform (on the large scale) intensity pro- 



file across the beam, though introducing into the beam 
small-scale, high-intensity variations or "speckles." In 
vacuum, speckles have characteristic size 8/^Ao (longi- 
tudinally) by /Ao (transversely), where / is the optical 
focal parameter and Aq is the laser wavelength. The 
scaling of SRS reflectivity Rsrs with laser intensity / 
in a solitary speckle in plasma has been measured [8] 
and found in the electron trapping regime /cAoe ^0-3 
(k is the EPW wave number and Aoe = ^/kB^^J^7^n^ 
is the Debye length for plasma of electron density ng 
and temperature Tg) to behave nonlinear ly, increasing 
sharply at a threshold intensity /th and saturating for 
/ > /th- The physics in this regime is governed by 
the growth of large-amplitude EPW that trap "resonant" 
electrons with speeds along the wave propagation direc- 
tion matching the wave's phase speed; this reduces local 
Landau damping [9], enhances instability growth, and 
lowers the EPW frequency [10 . At high intensity, trap- 
ping introduces variation in EPW phase velocity across 
the speckle and causes wave phase fronts to bend [TT1 - 
[14]. As EPW grow to large amplitude, secondary, non- 
linear processes have been proposed to break the phase 
fronts into small-transverse-scale filaments [fMHTB] that 
further contribute to nonlinear saturation. An effect of 
saturation, observed in simulations [5 , is the generation 
of hot electrons and back- and side-scattered light waves 
that propagate obliquely out of hot spots and enhance 
SRS growth in neighboring speckles through larger SRS 
seed levels and reduced EPW Landau damping. At high 
gain in two spatial dimensions (2D), this coupling enables 
networks of speckles to self-organize [5 and exhibit emer- 
gent behavior where reflectivity exceeds that of the sum 
of contributions from individual, non-interacting speck- 
les. The nonlinear nature of SRS in this regime is robust, 
with a threshold at modest laser intensity, > 10^^ W/cm^ 
for NIF laser conditions where /cAoe ^ 0.3 and the high- 
est levels of backscatter are found [17 . 

Our intent here is to show that the use of STUD 
pulses [2], effective for controlling LPI over long time 
scales in the fluid regime of instability evolution [2]-[4], 



may also inhibit EPW growth in the highly nonlinear, 
kinetic regime. STUD pulses deliver laser power in a se- 
quence of pulses or spikes on the instability growth or 
hot spot crossing time scale (sub-ps, typically, for SRS) 
with randomized laser speckle patterns in between one or 
more successive spikes. By introducing on-off sequences 
of pulses and by spatially scrambling locations of hot 
spots, reinforcing processes within a hot spot and the 
interconnectivity between hot spots that leads to large 
instability growth are disrupted. STUD pulses introduce 
degrees of freedom that can be optimized [2 . These in- 
clude the ratios Lhs • ^int * ^spike, where the inter- 
action length is Lint = 4Ln[a^/i4 + (^2/^2) ]^^^- Here, 
the density scale length is L^ — |Vlogne|~^, /14 is in- 
tensity in 10^^ W/cm^, V2 and a;2 are the local Landau 
damping rate and frequency of the EPW, respectively, 

a^ = 1.42A[r'"li:L'°° '""I 1^1/^2! (ne/ner)-\ and V^jV^ is 
the ratio of group velocities of scattered SRS light to 
the EPW [2 . Spike length Lgpike is the distance trav- 
eled by scattered light during the 'on' time Tgpike and 
^HS ^ 4/^Ao = 90 /im is the characteristic size of a 
hot spot in our plasma. Other degrees of freedom are 
duty cycle (the ratio of Tgpike to on+off time) , the spatial 
scrambling rate xriscram (how many identical spikes be- 
fore the RPP pattern changes) and the "jitter" (random 
variation of each Tgpike; all calculations in this Letter used 
0% jitter). Hence, "5000x1, 1:1:1/2" indicates a STUD 
pulse sequence with 50% duty cycle, 0% jitter, and a spike 
duration half as long as a hot spot crossing time in plasma 
where the time to cross I/jnt for the three-wave process 
is comparable to that of crossing the hot spot. Most of 
the results we present are for cases 5000x1, 1:1:1/2 and 
1:1:1 in strong to very strong nonlinear kinetic regimes 
(SRS gains of 6-11 at the average intensity). Note that 
configurations where 'on' time is much greater than 'off' 
time, e.g., 8000x1 or 9500x1, resemble the ISI model of 
beam smoothing [18 at the same bandwidth. 

To understand better the behavior of STUD pulses 
in the trapping regime, we have run VPIC particle-in- 
cell simulations [19 of a 2D plasma of size 500 x 80 /im 
in (x,z), with a laser beam polarized along y launched 
at X = as described in Ref. [6 . The laser has wave- 
length Ao = 0.351 /im and an RPP speckle pattern for 
/\8 speckles is used, approximating a NIF inner-cone 
beam in hohlraum plasma. The density has a gradient 
along X with ng = 0.12ncr at the center, changing from 
ibO.OlSncr to ±0.03ncr across the box, comparable to the 
Ln ~ mm encountered in NIF ignition hohlraums in re- 
gions of high SRS backscatter [17 . Taking V2 = 1^2^^^, 
as for Maxwellian plasma, in the /cAoe ^ 0-3 regime 
yields Lint ^ 95-99 /im for the range of intensities sim- 
ulated. We use 36864 x 4096 cells {Ax = 1.2ADe and 
Az = l.TAoe) and 256-512 electron macroparticles/cell; 
ions are stationary. The electrons have Tg = 2.6 keV 
(/cAoe = 0.3). The STUD pulse speckle patterns are 
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FIG. 1. (Color) Ey (left) and corresponding instantaneous back- 
scattered Poynting flux mdiyi (— Ey B z , 0) (right) over the 2D sim- 
ulation volume for three cases: (a) an RPP laser beam at aver- 
age laser intensity {/) = 5 x 10^^ W/cm^ (top), (b) a 5000x1, 
1:1:1/2 STUD pulse beam at time-averaged incident laser intensity 
{/) = 5 X 10^^ W/cm^ (center); and, (c) the same STUD pulse 
beam, but with {/> = 3.2 x 10^^ W/cm^ (bottom) (note logarith- 
mic scale on Poynting flux). The center panels are electrostatic fleld 
energy E^ for the leftmost 80 /im of each simulation, showing EPW 
wave amplitude correlated with instantaneous SRS backscattered 
Poynting flux. The inset is reflectivity vs. time for cases (a) (black) 
and (b) (blue); (c) (red) evinces negligible backscatter. The times 
shown are 1.6 ps (a) and 3.6 ps (b,c), chosen when large, bursts of 
SRS backscatter were present in (a) and (b). 



generated from pre-computed RPP phases for a wide 
beam, sampling 80-/im, non-overlapping segments for 
each STUD pulse. Each simulation employs the same 
sequence of speckle patterns to within an overall inten- 
sity modulation, allowing variation of intensity, duty cy- 
cle, and modulation period. (Statistical variation was as- 
sessed by altering the sequences of STUD pulses; ~10% 
relative Rsrs variation was found in a range of cases con- 
sidered.) The simulation boundaries absorb electromag- 
netic waves and reinject electrons as Maxwellian at initial 
temperature Tg. The simulations were run until apparent 
"steady-state" in time-averaged Rsrs^ 10-20 ps. 

Fig. 1 shows a comparison of three simulations: (a) 
(top row) is for an RPP beam with (/) = 5 x 
10^"^ W/cm^ (G = 11); (c) (bottom) is for a STUD 
pulse beam of time-averaged laser intensity (/) = 3.2 x 
10^"^ W/cm^ (G = 7.5). Linear SRS gains G are com- 
puted from G = 47r(7o/a;o)^ (27rLn/Ao)^(n/nc)""^(l - 



^1^2/70)7 where (70/^0^0) = 0.0043v7i4 A[f" '^ ui is the 
damping rate of the daughter light wave, and g{n) = 
^/Y^^2^^ [{1 / ^/n) — 1] [2. Accounting for backscat- 
ter loss, (a) and (c) have comparable net time- averaged 
power incident on the left boundary, though (c) has only 
64% of the incident time-averaged laser power. Case 
(b) (center) is for a STUD pulse beam at the same 
time-averaged incident laser intensity as (a): (/) = 5 x 
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FIG. 2. (Color) Hot electron flux per unit length <l> vs. elec- 
tron energy {Ek) for the three simulations in Fig. 1. Shown are 
trapped particle fluxes, obtained by subtracting contributions from 
a Maxwellian. Fluxes are measured on boundary regions, as in- 
dicated by the colors (c.f. the simulation box above, drawn to 
scale): z = ±40 /im, < a: < 250 //m (black); z = ±40 /im, 
250 < a: < 500 /im (red) and x = 500 /xm (blue). 




-30 -20-10 10 20 30 

Angle (deg.) 

FIG. 3. (Color) Angular distribution of the time-averaged 
backscattered light power for case (a) (black, dashed), (b) (dot- 
dashed, blue), and (c) (red) as in Figs. 1 and 2. The spectra for 
(b) and (c) evince lower backscattered light power, but at flnite an- 
gle with respect to the incident laser (cone angle |^| < 1/2/, shown 
by the vertical lines). 



10^"^ W/cm^ (G = 11). The leftmost panels show Ey (or 
the vacuum speckle pattern for the RPP case) . The right- 
most panels are instantaneous backscattered Poynting 
flux iRdix{—EyBz,0). Case (a) evinces continual bursts 
of self-organized backscatter with peak Rsrs > 1- In (c), 
no self-organization is seen in backscattered light or lon- 
gitudinal electric field. Case (b) is intermediate, showing 
quiescent periods of low backscatter punctuated by oc- 
casional episodes of partial self-organization when large- 
amplitude speckles (/ > 10 (/)) exhibit large-amplitude 
EPW and secondary processes, such as obliquely side- 
scattered light, occur at sufficient amplitude to seed SRS 
in otherwise stable regions of plasma (seen in the finite 
backscattered SRS Poynting fiux across the left of the 
box). The instantaneous Rsrs at the left boundary is 
shown in the inset for cases (a-c) (black, blue, and red 
curves, respectively); the times plotted are 1.6 ps for the 
RPP (during the first large SRS burst), and 3.6 ps for 
the STUD pulse simulations [during the first, large SRS 
burst in (b)]. The central panels are E'^ over the left- 
most 80 /im of the volume and indicate EPW amplitude 
correlated with the large bursts of SRS in (a) and (b). 

In Fig. 2, we compare for cases (a-c) time-integrated 
hot electron fiux per unit area exiting the simulation. 
The black curves are fiuxes leaving the ±z boundaries 
from the left half of the simulation volume, the red 
curves, leaving ±2: from the right of the volume, and 
the blue curves, leaving from the -\-z boundary. Prior 
work showed that large fiuxes of tail electrons leaving 
the left of the volume (i.e., large black curves) are sig- 
natures of large- amplitude EPW with ensuing nonlinear 
self- focusing and filament at ion and, ultimately, collective 
behavior among speckles [H [7] . The three cases evince 
elevated distribution function tails as a consequence of 



trapping, though the RPP case traps not only far more 
tail electrons [60 x more than (c), 6x more than (b)], but 
also shows far more side-scattered hot electrons exit near- 
est the laser entrance; moreover, hot electrons at very 
high energy {Ek > 100 keV) are present (absent for the 
STUD pulse beams). The use of STUD pulses has there- 
fore decreased the number of hot electrons exchanged 
laterally among laser speckles, a key part of inter-speckle 
self-organization [71 and a possible contributor to cap- 
sule preheat in ICF experiments. In Fig. 3, we compare 
the angular spread of SRS light for the same cases as 
above. The use of STUD pulses leads to a dramatic 
overall reduction in SRS power (and hence, amplitude 
of the SRS seed in neighboring speckles). As with the 
RPP, the angular spread is finite, with most of the power 
falling outside the incident laser cone \0\ < 1/2/ shown 
by the vertical lines. While the existence of coherent, 
oblique cones of backscattered light is not unique to this 
nonlinear regime — indeed, they appear in paraxial mod- 
els with diffraction [2|[20l — additional side-scatter results 
from trapping and EPW filamentation [S] [16] that is ab- 
sent in ffuid models; the use of STUD pulses reduces 
dramatically the levels of such side-scatter. 

Finally, in Fig. 4, we show the dependence of Rsrs 
on time-averaged incident laser intensity (left) and lin- 
ear gain at the average intensity (right) for RPP and 
STUD pulses over a range of plasma and laser condi- 
tions. The use of STUD pulses reduces dramatically 
Rsrs compared with RPP and ISI-like beams with the 
same time- averaged laser power. This is so even in cases 
of very high linear gain. As seen from comparison of 
the Rsrs from the ISI-like points (the 8000x1, 1:1:1/2 
and 9500x1, 1:1:1/2) and 5000x1, 1:1:1/2 cases, "heal- 
ing time" is key: it is not enough to simply add band- 
width and spatial scrambling. By optimizing this healing 



time for given "on+off" time and time- averaged power, 
STUD pulses may be optimized to significantly outper- 
form ISI. From comparison of the 5000x4, 1:1:1/2 and 
the 5000x1, 1:1:1/2 cases, we find that spatial scram- 
bling of the locations of the hot spots is also necessary to 
avoid effects of recurrence and correlation among suc- 
cessive hot spots [21]. Also, for fixed "on+off" time 
and time-averaged power, lengthening 'off' time requires 
shortening Tgpike and increasing the average speckle in- 
tensity correspondingly. Taken to an extreme, this can 
lead to enhanced trapping and associated EPW nonlin- 
earity, evidenced by the 2000x1 datum shown in Fig. 4 
[which also has significant hot electron sidescatter (not 
shown) compared with the 5000x1, 1:1:1/2 case at the 
same average power]. 

Examination of velocity distribution functions and 
EPW amplitudes shows strong trapping and only mod- 
est EPW damping between cycles. This trapping mod- 
ifies I/iNT and suggests that possible threshold behavior 
may arise when L^pike becomes less than Lint and when 
SRS goes from strong to weak damping. Consider the 
two 5000x1 and the 9500x1 STUD pulse cases at the 
highest intensity (G = 11). In the former, SRS in the 
largest amplitude hot spots (/ > 10 (/)) would be in the 
weak damping limit (WDL) if one applies the inferred 
U2 from simulations (~ 0.1i^2^^^), and Lint ^120 /im. 
The 1:1:1/2 STUD pulse case, with the lowest Rsrs^ has 
Lspike ~ 0.37LiNT for these maximal speckles, i.e., STUD 
pulses much shorter than Lint- In contrast, the 1:1:1 
case has Lgpike ^ Lint- In the 9500x1, 1:1:1/2 ISI-like 
case, while Lint = 0.5Lspike, reduction of U2 causes the 
SRS to go to the WDL for an average intensity speckle, 
with correspondingly large Rsrs- 

We have shown that SRS refiectivity may be lowered by 
more than an order of magnitude with the use of properly 
designed STUD pulses in settings where EPW trapping- 
induced nonlinearity is prevalent. This reduction stems 
from arresting large-amplitude EPW that give rise to 
cooperative behavior among laser speckles through the 
exchange of hot electrons and backscatter SRS waves, 
thus disallowing their self-organization. The substantial 
promise and generality of the STUD pulse technique [2] 
to a range of settings, including SRS in the strongly non- 
linear, trapping regime considered here, would seem to 
impel serious consideration for how STUD pulses might 
be achieved in future ICF laser systems such as the Green 
option on the NIF or next-generation high-repetition-rate 
laser systems. 
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FIG. 4. (Color) SRS reflectivity plotted as a function of average 
incident intensity (left) and linear gain (right) for RPP and STUD- 
pulse beams for a variety of laser and plasma conditions, showing 
the eflicacy of STUD pulses for reducing Rsrs- The black points 
(marked x and A) are for RPP beams for two different density 
variations across the box in x: zbO.OlSricr and zbO.OSricr, respec- 
tively. The red (EH) and blue (o) curves are for STUD pulses with 
5000x1, 1:1:1/2 and 5000x1, 1:1:1 (i.e., twice the "on" time) for 
the latter density proflle. The magenta datapoints labeled "ISI" are 
for STUD pulses with 8000x1, 1:1:1/2 (left) and 9500x1, 1:1:1/2 
(right) and indicate little advantage over RPP for the conditions 
shown. The green (A) is for STUD pulses with 5000x4, 1:1:1/2, 
i.e., scrambling every 4 pulses. The (+) is for 2000x1, 1:1:1/2. 
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